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ABSTRACT 



A series of reacting flow tests of a solid fuel ramjet 
were conducted with different air inlet swirl. Cold flow 
measurements of the air flow at the exit plane of the tube- 
in-hol^ injector were used to determine the Swirl Number for 
each configuration. Regression rates of an HTPB type fuel at 
low air mass flux were found to increase with swirl vane 
angle (and Swirl Number) up to 30 degrees. 
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I 



INTRODUCTION 



A. SOLID FUEL RAMJET 

Ramjet technology has taken three distinct approaches to 
meet the requirements of an ever-expanding mission responsi- 
bility “placed upon weapon systems of today. These three 
concepts are the Liquid Fuel Ramjet (LFRJ), Solid Fuel Ram- 
jet (SFRJ) and the Ducted Rocket (DR). Each system has 
distinct advantages over the other. The LFRJ has the 
capability to operate over widely different levels of thrust 
and efficiency through use of variable fuel control. The 
SFRJ is lighter, more compact, and is much less complex 
than the LFRJ due to the absence of fuel tanks, fuel 
delivery and control mechanisms. The DR retains some of 
the benefits of the SFRJ, but must carry oxidizing materials 
within the fuel to maintain combustion. 

Because of the simplicity with which the SFRJ operates, 
and its relative compactness, it is desireable to develop 
a simple thrust control for the solid fuel ramjet motor. To 
understand how this problem may be approached, it is first 
necessary to have a working knowledge of the basics of ramjet 
technology. 

As stated earlier, the ramjet is a very simple engine. 
Since there is no compressor, the ramjet uses supersonic 
velocities to compress the air flow. Figure 1 depicts the 
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Figure 1. Schematic of Simple SFRJ 
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INLET 



basic SFRJ. The three components of the basic SFRJ are the 
supersonic diffuser or inlet, the combustor, and the 
exhaust nozzle. 

As supersonic flows enter the SFRJ diffuser, several 
significant events take place. Standing shocks are present 
at the air inlet. An increase in static temperature and 
decreas?e in flow velocity occurs through the shock waves. 
Combustion of the solid fuel increases the mass flow slightly 
and the temperature is again increased. This larger mass 
flow at higher temperature is then expanded through the 
exhaust nozzle. A detailed description is given by United 
Technologies Corporation, Chemical Systems Division TRef. 1]. 
It is the increase (or decrease) of the fuel prolysis rate or 
regression rate (r) that determines the thrust and effects 
the efficiency of the SFRJ. 

There are many methods that can be used to affect the 
regression rate of a solid fuel. One is to change the fuel 
composition. Another is to change the environmental effects 
on the fuel, such as the chamber pressure fPc), inlet air 
temperature, and the mass flux through the port of the fuel 
grain. A typical regression rate formula is given by Boaz 
and Netzer [Ref. 2] 

for motors not using bypass air where 
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c = constant 

Pc = chamber pressure 

G = mass flow per unit area in the fuel port 
Ti = inlet air temperature 

The SFRJ is self-throttling through the effect of G on 
r. Hovfever, the self-throttling is not sufficient for many 
missions. Alternate means for varying r are therefore 
desired. Changing G can be accomplished by either changing 
m (such as in bypass designs) or by changing the flow area. 
Variable bypass is one possible means for controlling r. 

Ko [Ref. 31 investigated the effects of inlet air swirl 
on the regression rate of the solid fuel grain. Swirl 
imparts an angular component of mass flux into the flow field 
which can effect the residence time of the flow by increasing 
the effective length of the combusion section of the ramjet 
motor. It can also affect the convective heat transfer to 
the fuel surface and the mixing rates of fuel and oxygen. It 
was found that swirl may be an effective means for regression 
rate control, although only limited data were obtained. 

B. SWIRL NUMBER 

The introduction of swirl into a flow causes difficulty 
in characterization of the flow field. There is no longer 
only the axial component of the flow but now radial and 
tangential components as well. The Swirl Number (S) is one 
way to characterize the flow. Lilley and others fRef. 4] 
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have done considerable work with swirling flows and have set 
forth the non-dimensionalized swirl number (for specific 
assumptions about the inlet swirl behavior) as 



G (d/2) 

where ^ 

<5^ = tangential flux of swirl momentum 

Gx = axial flux of momentum 

d/2 = equivalent inlet radius 

This thesis investigated the effects of varying the inlet 
air flow angle on the combustion process of a non-bypass 
SFRJ. This was accomplished by using tube-in-hole injectors 
with vanes of 0, 15, 30, and 45 degrees. Cold flow 

measurements were made to determine inlet flow angles, and 
reacting flow tests were conducted to determine the effects 
on fuel regression rate. 
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II. DESCRIPTION OF APPARATUS 



A. RAMJET MOTOR 

The ramjet motor used in this investigation was identical 
with one used in earlier invesigations at he Naval 
Postgratduate School [Refs. 2,3]. Figure 2 is a schematic of 
the head-end assembly, step insert section, fuel grain and 
nozzle. 

The head-end assembly used a wedge to turn the primary 
inlet air from 90 degree plenum dumps. The igniter torch and 
and ignition gas were injected into the recirculation zone 
just aft of the step insert. 

The step insert section acted as the air inlet, which 
was modified with a tube-in-hole injector. The step height 
allowed for sudden expansion of the flow into the combustion 
section as depicted in Figure 3. This sudden expansion 
allows a flame to be stabilize behind the step by creating a 
recirculation zone. The swirl elements and tube were 
fabricated from stainless steel. Each swirl blade was a 
cambered airfoil attached at the hub with silver solder. The 
vane tips were ground for minimum clearance at the tube inner 
wall, and then silver soldered along the downstream side of 
the blades. Figure 4 shows the tubes with the swirlers used 
in this experiment. Tubes with 15, 30, and 45 deg. (from 
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. Depiction of Experimental SFRJ Motor Assembley 
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Figure 3. Step Inlet Configuration with Tube-in-hole 
Injector 
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Figure 4. 15, 30, and 45 Deqree Swirl 

Blades Mounted Inside Tubes 
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the axial centerline) swirl vane elements and a tube without 
a swirler were used. 

CSD-18818 (an HTPB based system) was the solid fuel used 
for the hot runs. The average length and internal diameter 
of the fuel grains used were 12.0 and 1.7 inches 
respectively. The fuel grain was placed between the head 
assembljj^^ and an aft mixing chamber. The latter contained the 
exhaust nozzle and was connected with steel rods which bolted 
to the head-end assembly. The entire motor was mounted on 
the thrust stand. Figure 5 is a photograph of the thrust 
stand a the test cell. The control panel was used to 
initiate the runs and the Hewlett-Packard 9836 Computer and 
3054A Data Acquisition/Control system sequenced events during 
the runs. 



19 




Figure 5. Experimental Thrust Stand and SFRJ 
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Figure 6. Schematic of Air and Gas Supply System 
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III. EXPERIMENTAL PROCEDURES 



A. CALIBRATION 

Calibration of pressure transducers was conducted 
at the beginning of each test day utilizing a dead-weight 
tester. Table 1 is a list of the transducers calibrated and 
the calibration pressures/load. 



TABLE 1 

TRANSDUCER CALIBRATION PRESSURES/LOAD 



Pa 


500 


psi 


Pc 


100 


ps i 


Ph 


100 


psi 


Phf 


400 


psi 


Pho 


500 


ps i 


Pif 


250 


psi 


F 


100 


Ibf 



B. DATA EXTRACTION 

A Honeywell 1508B Visicorder was used to obtain analog 
readings from pressure and thrust transducers. Chromel- 
alumel thermocouples were used for measurements of all gas 
temperatures. A Hewlett-Packard 9836 Computer and 3054A 
Automatic Data Acquisition/Control Stytem were utilized to 
extract digital information. 
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C. HOT FIRING 



Fuel grain length (L), average port diameter (Dp) , and 
initial weight (W) were measured for each run. The inlet was 
then mounted on the ramjet motor and preparations were made 
for the ensuing runs. 

Air and ignition gas pressures were calculated using the 
choked "flow formula (assuming Cd = 0.97) 




( R Tt) ^ 



Pressures were set using dome loaders and mass flow rates 
were set using the computer system as the controlling device. 
The computer was then readied for data acquisiton and each 
test was initiated from the control panel. Table 2 lists the 
inlet geometry and test condintions for each run. 

Upon completion of each run, the fuel grain was weighed 
and weight loss recorded. Actual motor burn time (tb), air 
time (ta), ignition time (ti), pre-ignition chamber pres- 
sure and thrust (Pa and Fa), and average chamber pressure 
and thrust (Pc and F) were obtained from the traces recorded 
on the Visicorder. The change in fuel grain weight is be 
expressed as 

AW = I (d| - d2)Lpg (4) 

Solving for the average final port diameter. 



d 



f 





(5) 
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TABLE 2 



* 



** 



INLET GEOMETRY AND TEST CONDITIONS 



Run 


Tube Length 


Swirl Angle 


ta, ti , tb. 




{ in. ) 


(deg. ) 


{ sec. )* 


1 


2.5 


0 


6,1, 6,4 


2 


2.5 


0 


6,1, 6,4 


3 


2.5 


15 


6,1, 6,4 


4 * 


2.5 


15 


6,1, 6,4 


5 


2.5 


30 


6,1, 6,4 


6 


2.5 


30 


6,1, 6,4 


7 


2.5 


45 


6,1, 6,4 


8 


2.5 


45 


6,1, 6,4 



Remarks ; 



Unsuccessful first two attempts. Cooled motor, 
reweighed grain, cleaned ignition gas choke and 
torch. Third try successful. 

pre-ignition airflow time, ignition time, burn time and 
purge time respectively 



24 



and substituting into 



r = 



(df - d. ) 

2 ^ 



( 6 ) 



yeilds the average fuel regression rate. 

Average values for m(air), m(02), m(C ) , 1^, and 
in(ign f.uel) were then calculated from the digital printout 
for the run. The fuel mass flow rate was calculated by 



m, 



A 

“V 



(7) 



The Naval Weapons Center (NWC) China Lake, Ca., 
Propellant Evaluation Program (PEPCODE) was used to obtain 
theoretical values for the adiabatic combustion temperature 
and combustion gas properties, Y and R, using Pc, Ti , 
heats of formation, and elemental composition by weight of 
the fuel grain and injected gases used during each run. Ko 
[Ref. 3] gives a complete explanation of the procedures used 
to obtain the values needed for data reduction. 



D. COLD FLOW PROCEDURES 

A hemispherical, three-hole stagnation pressure probe 
was used to measure flow velocities at the exit plane of the 
tube-in-hole injector. The probe was connected by tubing to 
a pressure transducer. The transducer in turn, was connected 
through a controller to a voltmeter where stagnation pres- 
sure readings in millivolts were displayed. A Scanivalve 
allowed for cycling through the three pressure readings to 
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find the angle at which the flow was eminating from the 
tube-in-hole injector exit plane. This was accomplished by 
finding the maximum value for the center port and nulling the 
left and right values. The velocities between the tube and 
the inlet wall were also measured. 

The probe could be rotated about its own axis. It was 

w 

mounted on a vernier translation stand to allow for movement 
across the flow in 0.05 in. increments. The probe itself 
could be translated on its own axis as well as rotated. 

This facilitated taking measurements in the y-z plane with 
rotation about the y-axis. Figure 7 is a photograph of the 
probe mounted on the thrust stand in this configuration. 

With the angle of flow taken at each grid point, the probe 
was then rotated from the y-axis to the z-axis. Figure 
8 depicts the probe head and the degrees of freedom in which 
it could be moved. 

The pressure readings in millivolts and the angles at 
which the the maximum values occurred were then input into a 
program that caluclated x, y, and z components of the flow 
velocity, and transformed them from the Cartesian coordinate 
system to the cylindrical coordinate system. This produced 
axial and tangential velocity components (assuming that the 
net radial component was equal to zero) . The computer 
program and examples of values calculated are given in 
Appendix A. 
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Figure 7. Photograph of Pressure Probe on Stand 
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Y-Z TRANSLATION 

Y-axis rotation 



Figure 8, Pressure Probe Head and Degrees of Freedom 
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Upon completion of characterization of the flow at each 
point in the grid, the amount of swirl in the flow was then 
characterized using Swirl Number (S) as previously stated 
in equation 2. The tangential flux of momentum is 



s = 



JL 

2 



u w 
o mo 



(d/2) 



(8) 



and th^ axial flux of momentum. is 

P _ JL P (1- (g/2)M 

“ 2 ° 

Making the same assumptions as Lilley [Ref. 4] of solid body 
rotation (a flat upstream axial velocity profile, and an 
increasing swirl velocity from 0 at the hub to w^ at the 
outer edge), the analysis yeilds the Swirl Number given as 



S = 



G/2 

(1 - (g/2) 2') 



( 10 ) 



where 

w 

( 11 ) 

mo 

is the ratio of the maximum tangential velocity to axial 
velocity at the exit plane of the tube-in-hole injector. 

In an attempt to more accurately characterize the swirl 
in the present apparatus, a mass averaged momentum ratio was 
calculated by taking the mass average tangential velocity 



w = 



Z m w 
Z m 



( 12 ) 
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and the mass average axial velocity is given by 



u 



Emu 
E m 



( 13 ) 



Dividing equation 12 by equation 13 yeilds S(avg). 
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IV. 



RESULTS AND DISCUSSION 



A. INITIAL TESTS 

A series of initial hot runs was completed to test the 
integrity of the tube-in-hole injector design. During this 
phase vt was found that initially designed tubes, which 
extended 0.25 inch beyond the aft end of the inlet, did not 
allow ignition of the motor. It was suspected that swirling 
flow was interferring with the recirculation zone, destroy- 
Ong the flame stabilization area. The tubes were lengthened 
so that the swirl dump was 1.5 inches beyond the inlet step, 
which solved the problem. This length was considerably less 
than the reattachment length of approximately 3.0 inches. 

Burn time for the motor was initially set for 10 seconds, 
which almost resulted in for burn-through to the wall. Burn 
time was therefore reduced to 8 seconds, which worked well 
for non-swirl test runs. Burn-through occurred at 8 secs, 
with 30 degrees of swirl, so the burn time was further re- 
duced to 6 seconds and the remaining hot runs were conducted. 

B. COLD FLOW ANALYSIS 

A summary of results of the cold flow studies is given in 
Table 3. There was an increase in the difference between the 
maximum and mass averaged axial velocities from 8% for no 
swirl, to 43% for 45 degrees of swirl. The same was true 
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for the difference between the maximum and mass averaged 
tangential velocities, with 21% difference at 15 degrees and 
43% difference at 45 degrees. This could have been a 
function of the increased blockage within the tube. The 
velocities between the tube and inlet wall increased as the 
vane angle increased, which was an indication that blockage 
was occurring. The Swirl number, as calculated by Equation 
10 (based on solid body rotation), increased with vane angle 
to 30 degrees and then decreased. The mass averaged value 
(Equation 12 divided by Equation 13) increased continually 
with vane angle. 

TABLE 3 

COLD FLOW RESULTS 



Blade angle 
(deg) 


0 


15 


30 


45 


u ( max ) 

( ft/sec) 


475 


616 


759 


820 


w( max ) 

( ft/sec) 


0 


309 


521 


525 


u ( avg) 

( ft/sec) 


434 


509 


560 


462 


u(outer flow) 
( ft/sec) 


384 


393 


503 


589 


w(avg) 

( ft/sec ) 


0 


243 


345 


296 


S ( avg ) 


0.0 


0.478 


0.616 


0.640 


S 


0.0 


0.335 


0.649 


0.544 
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C. FUEL GRAIN PHYSICAL CHARACTERISTICS 



The fuel grains used without swirl exhibited smooth 
burning surfaces with the port diameter slowly increasing to 
a point approximately 2.0 inches beyond the step. There- 
after, the port diameter was approximately constant. The 
grains used for the 15 degree swirl showed the same charact- 
eristic's. However, four areas that corresponded with the 
four vanes of the swirl element exhibited slightly indented 
burn regions extending from 2 inches beyond the step to 4 
inches, slowly blending into a concentric surface. The 30 
degree grains showed a marked increase in the indented burn 
regions, and they extended to 5 inches beyond the step. The 
45 degree grains were even more pronounced in the indented 
burned region but not as elongated. The angles of the 
indented burned regions corresponded to the respective swirl 
vane angles. Table 4 lists the physical characteristics of 
fuel grains pertaining to the tests. 

D. COMBUSTION PROPERTIES 

Measured quantities for each hot run are listed in Table 
5. A plot of the regression rate vs swirl vane angle is 
given in Figure 9 and vs S and S(avg) in Figure 10. 
Regression rate increased rapidly for small amounts (15 
degree vane angle) of swirl. Further increases in swirl vane 
angles (or Swirl Number) had much less effect on f. The 
Swirl Number (Equation 10) correllated reasonably well with 
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regression rate. S(avg) appeared to be in ouestion for the 
45 degree vane, perhaps due to the large amount of blockage. 

A typical value of G {mass flux) in an SFRJ is 0.6 
Ibm/in^-sec without bypass. Fifty percent bypass would then 
yeild a fuel port G of 0.3. With r'^^ G***, an approximately 
25% reduction in r could be achieved. A 30 degree swirl, in 
the present configuration, would yeild only about half of 
this change in r. However, it does appear that variable 
swirl can be an effectiv method of regression rate control in 
the SFRJ. 



TABLE 4 

FUEL GRAIN PHYSICAL CHARACTERISTICS 



Run # 


Vane 

Angle 


L 


Dp 




(deg) 


(in) 


(in) 


1 


0 


12.010 


1.750 


2 


0 


12.000 


1.749 


3 


15 


12.008 


1.747 


4 


15 


12.014 


1.759 


5 


30 


12.010 


1.754 


6 


30 


12.060 


1.740 


7 


45 


11.929 


1.750 


8 


45 


11.915 


1.746 



Di 


Dth 


Ath 


AW 


(in) 


(in) 


(sq. in) 


(gm) 


0.92 


0.943 


0.6483 


197 


0.92 


0.943 


0.6604 


200 


0.92 


0.943 


0.6554 


219 


0.92 


0.943 


0.6735 


222 


0.92 


0.943 


0.6597 


233 


0.92 


0.943 


0.6525 


230 


0.92 


0.943 


0.5721 


221 


0.92 


0.943 


0.6648 


220 
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1 

2 

3 

4 

5 

6 

7 

8 

1 

2 

3 

4 

5 

6 

7 



TABLE 5 



MEASURED COMBUSTION PROPERTIES 



ti 

( sec) 


tb 

( sec) 


ta 

( sec) 


Pa 

(psia) 


r 

( in/ sec) 


Pc 

(psia) 


Ti 
(® R) 


0.82 


7.27 


6.0 


50. 3 


0.0236 


100.2 


1103 


0.84 


7.26 


6.1 


49.9 


0.0240 


107.8 


1125 


0.86 


7.19 


6.0 


49.5 


0.0267 


105.7 


1125 


0.78 


7.20 


6.1 


50. 3 


0.0266 


107.8 


1119 


0.73 


7.21 


6.0 


51.1 


0.0277 


108.6 


1122 


0.80 


7.22 


6.1 


51.1 


0.0274 


109.3 


1119 


0.79 


7.18 


6.8 


50.7 


0.0263 


107.8 


1122 


0.87 


7.27 


6.0 


49.1 


0.0276 


103.8 


1130 



m m 

air fuel 

( Ibm/sec ) 

0.503 0.0598 

0.506 0.0608 

0.495 0.0672 

0.510 0.0680 

0.508 0.0713 

0.513 0.0702 

0.503 0.0602 

0.487 0.0667 



m 


m 


0^ 


CH 


2 


4 


0.0190 


0.0046 


0.0190 


0.0047 


0.0191 


0.0047 


0.0193 


0.0048 


0.0191 


0.0047 


0.0194 


0.0049 


0.0194 


0.0048 


0.0192 


0.0049 



m 


m 


if 


tot 


0.0033 


0.593 


0.0335 


0.593 


0.0034 


0.589 


0.0034 


0.604 


0.0034 


0.606 


0.0033 


0.611 


0.0034 


0.591 


0.0036 


0.581 
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REGRESSION RATE 

0.0230 0.0235 0.0240 0.0245 0.0250 0.0255 0.0260 0.0265 0.0270 0.0275 0.0280 



I 

t 



REGRESSION RATE VS VANE ANGLE 
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REGRESSION RATE 

0.0230 0.0235 0.0240 0.0245 0.0250 0.0255 0.0260 0.0265 0.0270 0.0275 0.0280 



REGRESSION RATE VS SWIRL NUMBER 
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V. CONCLUSIONS AND RECOMMENDATIONS 



It has been found that the regression rate increases 
significantly with small amounts of swirl in the flow of 
inlet air, but that larger amounts of swirl yield less 
effect.*^ The Swirl Number S (as given by Equation 10), at 
the low mass flux values of approximately 0.25 which were 
used in this investigation, correlated reasonably well with 
regression rate. 

It is recommended that the effects of Swirl Number on 
regression rate be investigated with increasing values of 
mass flux to 1.0 Ibm/in -sec and with metalized fuels. 
Improvement of the 45 degree swirl element design should be 
also attempted, as well as an investigation of the proper- 
ties of combustion with swirl angles between 30 and 45 
degrees and beyond. 
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A£££NQIX A 

PROGRAM FOR COLD FLOW ANALYSIS OF SWIRLING INLET AIRFLOWS 

SWIRL NUMBER PROGRAM 
LT. WILLIAM H. CAMPBELLt JR. 

NAVAL POSTGRADUATE SCHOOL, MONTEREY 

y 

THIS PROGRAM IS USED TO COMPUTE THE AVERAGE SWiRL NUMBER OF A 
SWIRLING FLOW OF A TUBE IN HOLE INJECTOR FCR A SOLID ROCKET 
RAMJET AND TRANSFORM THE CARTESIAN VELOCITY COMPONENTS INTO 
CYLINDRICAL VELOCITY COMPONENTS. 



DEFINITIONS 



TI ( li 


= 


INLET AIR TEMPERATURE 






SONVEL (I ) 




SONIC VELOCITY 






OENSd ) 


= 


AIR DENSITY 






Y( 1,K) 


= 


Y-CARTESIAN COORDINATE OF THE 


EXIT 


PLANE 


Z( J,l) 


= 


Z-CARTESIAN COORDINATE OF THE 


EXIT 


PLANE 


K1 ( J,K) 


= 


radius to point Y(N),Z(M) 






RY( I, J,K) 


= 


PRESSURE VECTOR IN MILLIVOLTS 


ABOUT 


THE Y-AXIS 


RZ( I, J,K) 


= 


pressure vector in MILLIVOLTS 


ABOUT 


THE Z-AXIS 


ALFA! I ,J,K) 


= 


ANGLE MEASURED OF RY(N,M), POSITIVE 


ABOVE THE X= 


BETA! I ,J,K) 




ANGLE MEASURED OF RZ(N,M), POSITIVE 


LEFT OF Y=0 


All I, J,K) 


= 


X COMPONENT OF RY(N,MI 






A2( ItJtK) 


= 


X COMPONENT OF RZ(N,M) 
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6(1, JfK) 
C(I,J,K) 
XCOMPi I ,J ,K)= 
RT GT(I,J,i^) = 

RXZ(I,J,K) = 

RXY(I,J,KJ = 

GAMMA ( I, J ,K)= 
OELTA( I,J ,K)= 
PHI(I,J,K) = 
TCOMP( I, J ,K) = 

RCOMP( I,J ,K)= 
PTX(I,J,K) = 
PTT(I,J,K) = 
U(I,J,K) 
^(I,J,K) 

SUMl 



SUM2 

SJM3 



Y COMPONENT OF RZ(N,M) 

Z COMPONENT OF RY(N,M) 

TOTAL X COMPONENT OF PRESSURE IN MILLIVOLTS 
TOTAL VECTOR OF PRESSURE IN MILLIVOLTS 
IN MILLIVOLTS 

TOTAL VECTOR COMPONENT OF PRESSURE IN THE X-Z 
PLANE IN MILLIVOLTS 

TOTAL VECTOR COMPONENT OF PRESSURE IN THE X-Y 
PLANE IN MILLI VOLTS 
ANGLE BETWEEN Y=0 AND RYZ 
ANGLE BETWEEN Y=0 AND R1 

ANGLE BETWEEN X=0 AND RXY, SWIRL ANGLE OF FLOW 
TANGENTIAL COMPONENT OF STAGNATION PRESSURE 
IN MILLIVOLTS 

RADIAL COMPONENT OF STAGNATION PRESSURE IN MILLIVOLTS 

X-COMPONENT OF STAGNATION PRESSURE IN PSI 

TANGENTIAL COMPONENT OF STAGNATION PRESSURE IN PSI 

AXIAL VELOCITY COMPONENT 

TANGENTIAL VELOCITY COMPONENT 

SUMATION OF THE GRID REPRESENTATIVE AREA 

MULTIPLIED BY THE G^ I 0 POINT TANGENTIAL 

VELOCITY W( I, J,K) 

SUMATION OF THE GRID REPRESENTATIVE AREAS 
SUMATION OF THE GRID REPRESENTATIVE AREA MULT. 

BY THE GRID POINT AXIAL VELOCITY U(I,J,KJ 
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UMAX< I ) 

WMAX( I ) 

AW( I ) 

AU( I) 

SAVd ) 
SMAXC I ) 
G( I ) 



THE MAXIMUM AXIAL VELOCITY FOR THAT RESPECTIVE 
SWIRL VANE ANGLE 

THE MAXIMUM TANGENTIAL VELOCITY FOR THAT 
RESPECTIVE SWIRL VANE ANGLE 
THE AVERAGE TANGENTIAL VELOCITY FOR THAT 
RESPECTIVE SWIRL VANE ANGLE 

THE AVERAGE AXIAL VELOCITY FOR THAT RESPECTIVE 

SWIRL VANE ANGLE 

THE RATIO OF THE AW( I ) TO AU(I) 

THE SWIRL NUMBER ACCORDING TO LILLEY ET AL 
THE RATIO OF THE WMAX ( I) TO UMAX! I) 



OELCARATION OF MATRICES FOR RESPECTIVE COMPONENTS OF THE PROGRAM 

REAL Y(l,8)» Z(8»l)» TI(4)» THETA(4)» SONVEL(^)» DENS ) » AW ( 4) ♦ 

+ RY(4»8t8), ALFA(4,8,8)t RZ(4,8»8), B ETA ( 4 ,8 t 8 ) , AO ( 4 ) , 

+ XCOMPI 4,8,8), Rl(8,8), RYZ(4,8,8), RTOT ( 4, 8 , 8 ) , SUMl ( 4 ) , 

+ PHI(4,8,8), RCCMP(4,8,8) , PTX(4,8,3), U(4,8,8), SUM2I4) 

REAL PTT(4,8,8), W(4,3,8), G(4), AREA(8,3), SUM3(4J, SUM4I4), 

+ Al(4,8,8), A2(4,8,8), B(4,8,8), C(4,8,8),RXZ(4,8,8), 

+ RXY(4,8,8), 0ELTA{4,8,8) , GAMMA { 4 , 8 , 8 ) , TCGMP ( 4 , 8 , 8 ) , 

+ WMAXI4), UMAXI4), SAV(4), SMAXI4) 

INCREMENTATION OF THE Y AND Z AXIS 
DO I N=l,8 
REA0I5, 10) Yd ,N) 
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1 



Z(N, ll»Y(l,N) 

CONT INUE 
DO 81 J = lt8 

R£A0(5,400) (AREA(JfK), K=l,8) 

81 CONTINUE 

COMPUTATION OF SONIC VELOC I TY , DENS ITY , AND INPUT OF VANE ANGLE. 
DO 2 1=1,4 

REA0(5,20) Tim, THETA! I) 

SONVEL ( n=SQRT{ 1 .4*32 .2*53 . 3*T I ( I ) ) 

DENS! I )=14.7/(53.3*TI(in 
2 CONTINUE 

THIS PORTION OF THE PROGRAM INPUTS AND OUTPUTS DATA RECORDED 
DURING THE COLD FLOW ANALYSIS OF DIFFERENT SWIRLER VANE 
GEOMETRIES.TH VALUE (8,8) FOR EACH ANGLE IS THE COMPONENT 
MEASURED BETWEEN THE INLET WALL AND THE TUBE-IN-HOLE INJECTOR 
OUTER WALL. 



DO 4 1=1,4 
SUMl m =0 
SUM2 (I) =0 
SUM3(I) =0 
SUM4 !I ) =0 
UMAX!! J =0 
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WMAX(I) =0 
DO 5 K = l,8 
DO 6 J=l,8 

READ(5,30) RY(I,J,K), ALFA(I,J,K), RZ(I,J,K), BETA(I,J,K) 
Al( I ,J,K)=RY( I, J,K)*CDS (ALFAd ,J,K) *3. 141/180) 

A2( I ,J,K)=RZ( I, J,K)»COS (BETAd , J,K)*3. 141/180) 

Bd , J,K) =RZ< I, J,K)*SIN (BETAd , J,K)*3.141/180) 

Cd ,J'»K)=RY(I ,J,K>«SIN(ALFA(I, J,K)*3.141/180) 

XCOMPd ,J,K) = Ald, J,K) + A2Ut Jf K) 

Rl( J,K)=SQRT( Y(1,K)*<=2 + ZU,1)*«2) 

RYZ d,J,K)=SCRT(Cd, J,K)**2+B( I,J,K )**2) 

RXZ ( I,J,K)=SQRT(Cd, J,K)**2+XCOMPd,J,K)=«‘*2) 

RXY < I,J,K)=SQRT ( XCOMP (I t J t K ) ♦* 2+B (I , J , K ) «*2 ) 

RTOKI, J,K) =SQRT(XC0MPd,J,K)««2+6( I , J , K ) «*2+C d , J , K ) **2 ) 
IF( ( Y(1 ,K) .EQ.O) .AND. (Z( J»l) .EQ.O) ) THEN 
0£LTA( I ,J ,K) =0 

ELSE IF ( (Z( J,l) .GE.0.05 ) .AND. ( Y ( 1 ,K) .EQ.O ) ) THEN 

DELTAdtJ»K)=90 

ELSE 

DELTA( I,J,K)=(ATAN(Z(J,l)/Yd,K) )) *180/3. 14 
ENO IF 

IF( Bd, J,K) .EQ.O) THEN 
GAMMA! I , J,K) =0 
ELSE 

GAMMA! I, J,K) =!ATAN!C! I , J , K )/ B ! I , J , K) ) ) *18 0/ 3. 14 
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END IF 

IF (XCOMPI I,J,K).EQ.O) THEN 
PHI(IfJfK)*0 
ELSE 

PH I(I,J,K) = ( ATAN(B( I, J,K)/XCOMP( I, J,K) ) )*180/3.1-V 
END IF 

TCONP( I, J,K ) = RY2(I, J,K)*COS( (93- (DELTA! I , J, K) +GAMMA( I , J,K) ) ) 
+ ♦3.141/180) 

RCOMP(I,J,K)=RYZ(I, J,K)^SIN( (90 -(DELTA! I , J, K ) +GAMMA ( I , J ,K) ) ) 
+♦3. 141/ 180) 

PTX ( I,J,K)=XCOMP( I ,J,K )/.0155 

U(It J,K)=SQRT( 2402. 7^TI (I )^(5^( ( (PTX ( I , J , <) / 14. 7) +1 ) ♦♦ .2857 1 
+ 1 ) ) ) 

PTT( I,J,K)=TCCMP( I ,J,K )/.0l55 

W(I ,J,K)=SQRT(2402.7^TI ( I ) ♦ ( 5* ( ( ( PTT ( I , J ,I0 / 14.7 ) + 1 ) ♦♦ 
+-!))) 

SUM1(I)=SUM1( I)+(AREA( J,K)^W(I ,J,K) ) 

SU'12(I)=SJM2( I)+AREA(J,K) 

SUM3(I )=SUM3( I)+(AREA( J,K)^U(I ,J,K) ) 

IF (U(I,J,K).GT.UNAX(I ) ) THEN 
UMAX! I )=U( If JfK) 

ELSE 

UMAX! I )=UMAX (I ) 

END IF 

IF ( W( I , J, K).GT.WMAX (I ) ) THEN 



.28571 
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WHAX(I)=M( It JfK) 

ELSE 

WHAX( I )=WMAX (I ) 

END IF 
6 CONTINUE 

5 CONTINUE 

AW( IJ = SUM1( I)/SUM2( I) 

AU( I)=’5UM3( I)/SUM2(I ) 

SAV(I)=AW(I)/AU( I) 

G( I ) = WMAX( I )/UMAX( II 
SNAXd )=(G( I) /2)/(l-G(I )»*2) 

^ CONTINUE 

WRITE(6,410) (UMAXm, 1=1,4) 
WRITE(6,41l) (WMAX(I), 1=1,4) 
rtRITE(6,412) (AW(I), 1=1,4) 
WRITE(6,413) (AU(I), 1 = 1,4) 
WRITE(6,414) (SAV(I), 1=1,4) 
WkITE(6,415) (SMAX(I), 1=1,4) 

00 8 1=1,4 
i»RlTEI6,40) THETA(I) 

WRITE(6,50I (Y(1,L), L=l,3) 

WR ITE( 6,55 ) 

00 9 J=l,3 

WRITE(6,60) Z(J,l), (U(I , J,K) ,K=1,3) 
9 CONTINUE 
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3 



CONTINUE 
DO 12 1=1*4 
WRITE! 6*70) THETA! I) 

WRITE!6,50) ( Y( I , U » L = 1 , 8 i 
WRITE! 6,55) 

DO 11 J=l,8 

WRITE!6,60) Z!J,1), ! W ! I , J *K ) *K= 1 , 8 ) 

11 CONTfNUE 

12 CONTINUE 

DO 14 1=1,4 

WRITE!6,80) THETA!I) 

WRITE!6,50) 1 Y! 1, L) , L= 1, 8 ) 

WRITE ! 6,55) 

DO 13 J=l,8 

WRITE!6,90) Z!J,1J, !RYII,J,K), K=l,8) 

13 CONTINUE 
lA CONTINUE 

00 16 1=1,4 
WRITE!6,80) THETAII) 

WRITE!6,50) ! Y! 1 , L ) ,L= 1 , 8 ) 

WRITE! 6,55) 

00 15 J=l,8 

WRITE!6,120) Z!J,1), !PHI!I,J,K), K=l,8) 

15 CONTINUE 

16 CONTINUE 
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17 

18 



21 

22 



23 

24 



00 18 I-lt4 

WRITE(6#100) THETA(I) 

WRITE(6«50) (Y( 1*U »L=1»8) 

WRITE (6,55) 

00 17 J=l,8 

WRITE(6t90) Z(J«D* (RZ(L«J»K)» K=l,8) 
CONTINUE 
CONT I^IUE 
DO 22 I = 1,4 

WRITE(6,110) THETA(I) 

WRITE(6,50I (Y( 1,U ,L=1,8) 

WRITE (6,55) 

00 21 J=l,8 

WRITE(6,120) Z(J,1), (ALFA( I,J, KJ, K=l,8) 
CONTINUE 
CONTINUE 
00 24 1=1,4 

WRITE(6,130) THETA(I) 

WRITE(6,50) (Y(1,L), L=l,8) 

WRITE (6,55) 

00 23 J=l,8 

WRITE(6,120) Z(J,1), (3ETA( I ,J, K), K=l,3) 
CONTINUE 
CONT INUE 
DC 32 1=1,4 
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HRITE(6,140) THETA(I) 

WRITE(6,50) (Y(1,L), L=l,8J 
WRITE (6t55) 

DO 31 J=lf8 

WRITE(6,150) ZU,1)» (C(I,J,K», K=l,8) 

31 CONTINUE 

32 CONTINUE 

DO 34^1=1,4 

WRITE(6,160) THETA(I) 

WRITE(6,50) (Y<1,LI, L=l,8) 

WRITE! 6»55) 

DO 33 J=li8 

WRIT£(6,150) ZiJtUf (XCOMP ( I , J ,K) , K=l,8) 

33 CONTINUE 

34 CONTINUE 

DO 36 1=1,4 

WRITE«6,170) THETA(I) 

WRITE(6,50) ( Y( 1,LJ ,L=1,8) 

WRITE (6,55) 

DO 35 J=l,8 

WRITE(6,150) Z(J,l), (RTOT(I,J,K), K=l,8) 

35 CONTINUE 

36 CONTINUE 

DO 38 1=1,4 

WRITE(6,180) THETA(I) 
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WRITEC6f50J (Y(l,L), L = l,3) 

WRITE ( 6,55) 

00 37 J=l,8 

WRITE(6,120) Z(J,1), (PHI(I,J,K), K=l,3) 

37 CONTINUE 

38 CONTINUE 

00 42 1=1,4 

WRITE(6,190) THETA(I) 

WkITE(6,50) (Y(1,L), L=l,8) 

WRITE (6,55) 

00 41 J=l,8 

WRITE(6,150) Z(J,l), (A1(I,J,K), K=L,8) 

41 CONTINUE 

42 CONTINUE 

00 44 1=1,4 

WRITE(6,200) THETA(I) 

WRIT£(6,50) (Y(1,L), L=l,8) 

WRITE(6,55) 

00 43 J=l,8 

WRITE(6,150) Z(J,1), (A2(I,J,K), K=l,8) 

43 CONTINUE 

44 CONTINUE 

00 46 1=1,4 
WRITE(6,210) THETA(I) 

WRITE(6,50) (Y(1,U, L=l,8) 
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MRITE (6,55) 

00 45 J=l,8 

WRITE(6tl50) K=l,3) 

45 CONTINUE 

46 CONTINUE 

00 48 1=1,4 

WRITE (6,220) THETAt I ) 

WRITE(6,50) (Y( 1,L) , L = 1 ,8) 

WRITE (6,55) 

00 47 J=l,8 

WRITE(6,150) Z(J,1), (RyZn,J,K), K = l,8) 

47 CONTINUE 

48 CONTINUE 

00 52 1=1,4 

WRITE(6,230) THETA(I) 

WRITE(6,50) (Y(1,U, L = l,8) 

WRITE (6,55) 

00 51 J=l,8 

WRITE(6,150) Z(J,1), (RXZ<I,J,K), K = l,8) 

51 CONTINUE 

52 CONTINUE 

03 54 1=1,4 

WRITE(6,240) THETAl I ) 

WRITE(6,50) (Y( 1,L) , L = l,8) 

WRITE (6,55) 
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00 53 J=l,8 

WRITE(6,150) (RXY(I,J»K), K=l,8) 

53 CONTINUE 

54 CONTINUE 

00 57 1=1,4 

WRITE(6,250) THETA(I) 

WRITE (6,50) ( Y( 1,U , L=l ,8 ) 

WR ITE16,55) 

00 56 J=l,8 

WRITE(6,120) Z(J,l), (0£LTA( I,J ,K) , K=l,8J 

56 CONTINUE 

57 CONTINUE 

DO 59 1=1,4 

WRITE(6,260) THETA(I) 

WRITE(6,50J (Y(l,L), L = l,8) 

WRITE ( 6,55) 

00 58 J=l,8 

WRITE(6,120) Z(J,1), (GAMNA( I ,J ,K) , K=l,8) 

58 CONTINUE 

59 CONTINUE 
DO 62 1=1,4 

WRITE( 6,270) THETAd ) 

WRITE(6,50) (Y(l,L), L = l,8) 

WRITE ( 6,55) 

00 6 1 J=l,8 
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WRITE(6,150) Z(J,1), (TCOMP( I,J »K) , K=l,3) 

61 CONTINUE 

62 CONTINUE 

DO 64 I =lt4 

WRITei6,280) THETAd) 

WRITE(6»50) (YU,L)» L = l,8) 

WRITE (6,55) 

DO^'6 3 J = l,8 

WRITE(6,150) Z(J,l), (RCOMP( I,J ,K) , K=l,3) 

63 CONTINUE 

64 CONTINUE 

DO 66 I =1,4 

WRITE( 6,290) THETAd) 

WRITE(6,50) (Y(1,L), L=l,8) 

WRITE (6,55) 

DO 65 J=l,8 

WRITE(6,60) Z(J,1), (PTX(I,J,K), K=l,8) 

65 CONTINUE 

66 CONTINUE 

DO 68 1=1,4 

WRITE(6,310) THETAd) 

WRITE(6,50) (Y(1,L), L=l,8) 

WRITE (6,55) 

DO 67 J=l,8 

WRITE(6,60) Z(J,l), (PTT(I,J,K), K=l,3) 
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67 CONTINUE 

68 CONTINUE 

MRITE ( 6.330) 

WRITE(6,50) (Y(l,U, L = l,3) 

WRITE(6,55) 

00 73 J=l,8 

WRITE(6,60) Z(J,1J, (RKJ.KJ, K = l,8) 

73 continue 

10 FORMAT! F6. 3) 

20 FORMAT! 2F6. 3) 

30 FORMAT! 4F10.3) 

40 FORMAT! • 1* . *U! I ,J.K) FOR SWIRL VANE BLADE ANGLE F5. 1 , ' 

+ .'/ ) 

50 F0RMAT!9X,'Y* .8F3.3) 

55 FORMAT! '0* ,3X, *Z' ) 

60 FORMAT! F6,3,4X,8F 3.2) 

70 FORMAT! • 1* » *W! I ,J ,K) FOR SWIRL VANE BLADE ANGL E‘ , F5 . 1 , • 
♦ . V J 

80 FORMAT! • 1* , *RY! I, J.K) FOR SWIRL VANE BLADE ANGLE*, F5.1, 

+'.*/) 

90 FORMAT! F6. 3, 4X.8F 8.3) 

100 FORMAT! * 1*. *RZ! I, J.K) FDR SWIRLVANE BLADE ANGLE * ,F5. 1 , • 

+ . •/ ) 

110 FORMAT! * 1*. 'ALFA! I. J.K) FOR SWIRL VANE BLADE ANGLE*, F5. 

+ 1 ,* . •/) 
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120 

130 

140 

150 

160 

170 

180 

190 

200 

210 

220 

230 

240 

250 



FORMAK F6.3,4X,8F8.0) 

FORMAK *l*t*8ETAt I, J»K) FOR SWIRL VANE BLADE ANGLE*, F5. 

+ 1 ,* . •/) 

FORMAT! • 1*, 'C< I, J,K) FOR SWIRL VANE BLADE ANGLE • , F 5. 1 , * 
+ . •/) 

FORMAT! F6.3,4X,8F8.3) 

F0RMAT!v/1*,*XC0MP!I,J,K) FOR SWIRL VANE BLADE ANGLE*, F5 
+ .1, * .*/ i 

FORMAT! * 1*, *RTOT( I, J,K) FOR SWIRL VANE BLADE ANGLE’, F5. 

+ 1 ,* . */) 

FORMAT! * 1*,*PHI( I, J,K) FOR SWIRL VANE BLADE ANGLE*, F5.1 

FORMAT! * 1*, *A1( I, J,K) FOR SWIRL VANE BLADE ANGLE*, F5.1, 
*-*.*/) 

FORMAT! *1*, *A2!I,J,KJ FOR SWIRL VANE BLADE ANGLE*, F5.1, 

+*.*/) 

FORMAT! * 1», *B! I ,J,K) FOR SWIRL VANE BLADE ANGLE * , F 5. 1 , * 
+ . */ ) 

FORM AT! * 1* , *RYZ! I , J,K) FORSWIRL VANE BLADE ANGLE*, F5.1, 

♦*.*/) 

FORMAT! * 1* , *RXZ! I ,J,K) FOR SWIRL VANE BLADE ANGLE*, F5.1 

+,*.*/) 

FORMAT! * 1* ,*RXY! I ,J,K) FOR SWIRL VANE BLADE ANGLE', F5.1 

+,*.*/) 

FORMAT! * 1*,*DELTA!I,J,K) FOR SWIRL VANE BlADE ANGLE*, F5 
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+ . 1 , • .'/ ) 

260 FORMAK • 1* f'GAMMAd ,J,K) FOR SWIRL VANE BLADE ANGLE*, F5 

+ . 1 , • .'/ ) 

270 FORMAT! • 1* , 'TCOMP (I ,J,K) FORSWIRL VANE BLADE A NGLE * , F 5 . 



280 

290 

310 



330 

400 

410 

411 

412 

413 

414 

415 



FORMAT! • 1', *RCOMP(I,J,K) FORSWiRL VANE BLADE ANGLE', F5. 

+ 1 ,* . •/) 

FORMAtV' 1' , *PTX! I ,J,K) FOR SWIRL VANE BLADE ANGLE*, F5.1 

+ , '. * /) 

FORMAT! ' 1* , *PTT! I ,J,K) FOR SWIRL VANE BLADE ANGLE', F5.1 
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